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1 Introduction 



In a famous 1967 paper [0 Steven Weinberg asked the question: "What relations 
are imposed by current algebra upon the spectra of the 1 + and 1~ mesons?" Under 
certain conditions the answer was two sum rules involving the vector and axial- 
vector spectral densities. They are known as the Weinberg sum rules: 

I / — [Pv(s) — Pa(s)] = F% , (1) 
Jo s 

II / ds[p v (s)-p A (s)} =0. (2) 
Jo 

Assuming vector meson dominance and the KSFR relation these sum rules lead 
to the prediction that the p and a\ masses are related by m ai = y2m p , which is 
approximately valid. In this paper we ask two questions. The first one is: Given 
that QCD is now known to be the theory of the strong interactions, what extra 
information can we get from sum rules of the Weinberg type? 

The last fifteen years has seen a great deal of activity surrounding QCD at 
finite temperature. Of particular interest are the issues of deconfinement and chiral 
symmetry restoration at temperatures of the order of 160 MeV. Therefore, we are 
led to consider a second question: What are the implications of the approach to 
chiral symmetry restoration at finite temperature for sum rules of the Weinberg 
type? 

The status of the original Weinberg sum rules in the context of QCD sum rules 
was discussed by Shifman, Vainshtein and Zakharov [g] and then by Narison 
while a more up-to-date phenomenological analysis was performed by Peccei and 
Sola 0]. The two sum rules derived by Weinberg are very general, as he showed, 
and do not depend on specific details of the QCD Lagrangian. Higher order sum 
rules (involving more powers of s in the integrand) do depend on dynamics of chiral 
symmetry breaking in the vacuum. In section 2 we derive a third sum rule of the 
type of eqs. (1-2). This new sum rule involves the vacuum expectation value (VEV) 



2 



of a certain local four-quark operator. It can be obtained from the sum rule if we 
know the vector and axial-vector spectral densities accurately enough from exper- 
iment. It can also be obtained from lattice QCD; the chirality-violating structure 
of the operator helps here because its VEV has no short distance perturbative con- 
tribution. We perform a detailed analysis of all sum rules in section 3. We will see 
that they are restrictive enough to fill in gaps in the experimental data, allowing us 
to determine the spectral densities with quite some accuracy. 

There has been a lot of discussion in the literature and at conferences about 
the temperature dependence of hadron masses. Some calculations yield increas- 
ing masses, some yield decreasing masses, and still others yield masses that ei- 
ther increase or decrease depending on the quantum numbers of the hadron; see 
H 1) @) !> Hi an d the review |^0| . Clearly all these calculations are only approx- 
imate. In addition, the very notion of a mass at finite temperature must be very 
clearly defined, such as the screening mass or the pole mass corresponding to col- 
lective excitations. 

A common denominator of all studies of this type is the temperature depen- 
dence of correlation functions. It would be good if some general statements about 
these correlation functions could be made which rely on the fundamental properties 
of QCD at finite temperature. This is the aim in section 4. We generalize the 
original Weinberg sum rules, and the new one, to finite temperature. The first one 
(eq. (1)) generalizes to a sum rule involving only the longitudinal spectral density 
and depends on three- momentum. The second one (eq. (2)) generalizes to two 
separate sum rules, one involving the longitudinal spectral density and the other 
involving the transverse spectral density, both depending on momentum. At zero 
three-momentum they collapse to the same expression. In the vacuum there is no 
dependence on momentum because of Lorentz invariance, but at finite temperature 
there is a preferred rest frame, hence a dependence on momentum and on polariza- 
tion. We would like to point out here that probably the first discussion of Weinberg 
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sum rules at finite temperature was given by Bochkarev and Shaposhnikov 0] in 
the context of QCD sum rules and for zero momentum. 

The finite temperature sum rules can be used to constrain models or approxi- 
mations to QCD, and can help us to understand the approach to chiral symmetry 
restoration. Various possibilities will be considered in section 5. We should refer- 
ence here the early paper on phenomenology of the chiral phase transition in heavy 
ion collisions by Pisarski |11]. For recent discussion of the topic one can see [O]. 

We remark that throughout this paper we assume that the up and down quark 
masses are identically zero so that chiral symmetry is exact. Consideration of the 
impact of nonzero quark masses on the original Weinberg sum rules within pertur- 
bative QCD was done by Floratos, Narison and de Rafael [13]. 



2 Derivation of Zero Temperature Sum Rules from QCD 

We define the vector and axial-vector currents, 



ya 



q^(r a /2)q , 
Al = W5(T a /%, 



(3) 
(4) 



where r a /2 is the isospin generator. With this normalization the current algebra of 
charges obeys the equal time commutation relations 



Qvi Qv 


= ie abc Q c v , 


(5) 


Qvi Qa 


= ie abc Q A , 


(6) 


Q% Qa 


= ie abc Q c v . 


(7) 



We define the vector and axial-vector spectral densities in the usual way. They are 
positive definite quantities defined for positive s. 

< 0\V£(x)V?(0)\0 > = -7^3 / (sT ~ Pv(s) , (8) 



<OK(x)^(0)|0> 



(2tt) e 



<fp6{p u )e 



,0\ Jp-x 



p2 



p A (s) + F*8(s)p^ 



0) 



The dimension of the spectral densities is energy-squared. Note the explicit contri- 
bution of the pion to the axial- vector correlator. 

In this paper we work in imaginary time so that all distances are space-like, or 



Euclidean: x = t 



-t 2 . In this domain the spectral representation of the 



correlation functions looks as follows 



AD*(r) =<0\T[v ai "{x)V^(0) - A a ^(x)A b tM (0) 

§ ab 



\0> 



Att 2 t Jo 



dsy/s 3pv(s) - 3pa(s) - s F n 5(s) K\{y/sr) , 



(10) 



and 



AUS(t) ^<0\T\v a °(x)V b °(0) - A° a (x)A° b (0) 



lo > 



>ab 



4vr 2 



T JO 



ds 



;n) 



Notice that the integrands essentially involve the standard Feynman propagator for 
a particle of mass m which, in the Euclidean domain, is 



D(rn, 7~)f ree scalar 



m 

4-7T 2 '/ 



■ifi(mr) . 



(12) 



Exponential decay of the Bessel function K\ at large argument ensures convergence 
of such integrals for any QCD correlation functions, except probably at r = 0. In 
this sense, there is no difference between the Euclidean time representation |pX| ] 
and the Borel-transformed sum rules Q, in which the propagator is replaced with 
exp(— s/M 2 ), with Borel parameter M replacing Euclidean time r. 

The coordinate representation is more transparent and accessible to numerical 
methods, such as lattice calculations. Recent studies based on the instanton liquid 



model [14] and lattice QCD [15] have reported on the calculation ofaset ofEuclidean 



correlation functions, including vector and axial ones. Unfortunately, none of them 
has focused on their difference with sufficiently high accuracy, and therefore they 
are not discussed in the present work. 

We now come to the central idea behind the derivation of the sum rules: each 
sum rule corresponds to a particular term in the small-distance asymptotic expan- 
sion of the correlation function. 

In the limit t — > the product of currents can be expanded according to the 
operator product expansion (OPE), a very powerful means for connecting VEV's 
of quark and gluon operators to experimentally observable hadronic properties. 
The first terms in this expansion were first computed in Q. For the contracted 
polarization tensor the result is 



where, in the argument of the logarithm, \x « 1/r is the renormalization scale, 
and O p is a complicated four-quark operator. There is a similar expression for the 
correlator of two axial- vector currents but with a different four-quark operator O ai . 
For our purposes we only need their difference, which is given below. 

Since chiral symmetry breaking is a long wavelength phenomenon, at very short 
distances, or at very high energies, the difference between vector and axial-vector 
correlators should go to zero. Indeed, taking this difference one finds that all terms 
except for the four-quark operators in eq. (13) drop out. 

One can now look for consequences of this statement for the spectral density. 
Expanding the Bessel function in eq. (10) for small values of r we get 



35 ab 



1 + 



D a ^{r) = < 0|T V a »(x)V*(0) |0 > 
(r) <0| W 2 ^ 4 *V» 




hxfjiT) < 0\O p \Q > +■■■ (13) 




(14) 
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where C is Euler's constant. The OPE has no power divergence in r in the difference 
AD^*. Therefore the coefficient of 1/t 2 in eq. (14) must vanish. This is just the 
second Weinberg sum rule (eq. (2)). In the OPE framework it simply follows from 
the observation that the first covariant operators which are not chirality blind are 
four-quark ones which have dimension 6 or more. Similarly expanding eq. (11) for 
small r and applying the observation of chirality blindness we get 



00 ds r 

Pv{s) - PA (s) - sF%6(s) 



1 | s 

^4 + 4^2 



0. (15) 

/,! | 7- IT- | 

The first and second terms in the last square brackets reproduce the first and second 
Weinberg sum rules, respectively. 

The next term in the small r expansion is logarithmic. In eq. (14) we multi- 
ply the argument of the logarithm by p/p which we must do to match the OPE. 
Equating the coefficients of ln(/ir) in AD^(r) we obtain the third sum rule, 

III / dss [ Pv (s) - p A (s)] = -2tt < 0|a s O£|0 > . (16) 

Jo 

Here 

= {u Ll H a u L - d Ll H a d L ) {u R Yt a u R - d R Yt a d R ) , (17) 

where t a are the color SU(3) matrices and R, L stand for right and left-handed 
quarks. Note the appearance of the renormalization scale [i on the right side of 
this sum rule. Since the other side of the equation is expressed in terms of physical 
observables, it must be that a s (fi) times the four-quark operator is a renormalization 
group invariant. 

The numerical value of the VEV of this operator is unknown. The estimate 
suggested in [g] is based on the so called "vacuum dominance" hypothesis, which 
leads to 

< 0|O£|0 > = - < 0\uu\0 > 2 . (18) 
The accuracy of this estimate should of course be questioned, and various mod- 



els of chiral symmetry breaking [16] and lattice numerical calculations can be 



used for that purpose. Let us only add a comment on //-dependence here. If 
the vacuum dominance hypothesis is correct, then the VEV should be propor- 
tional to [ln(/i/Agc*£j)] 8 / 6 , the anomalous dimension of the quark condensate. (Here 
b=^-N c — |iVy comes from the Gell-Mann-Low function.) Since the power is close 
to 1, after being multiplied by a s (fi) ~ 1/ \n{p/ Kqcd) the right side of the third 
sum rule is nearly //-independent. Thus, at least concerning the //-dependence, this 
approximation can approximately hold. 

The regular (r-independent) term was not considered in the QCD sum rule 
context; it was first discussed in connection with point-to-point correlators in the 
coordinate representation by one of us . It is interesting to express it in terms of 
an integral over the difference in spectral densities, and it may be useful for lattice 
calculations. Dropping terms which vanish in the limit, we find 

3 f°° 

AZ>£(r-»0) = -lnOtr^jf ds s [p v (s) - p A (s)} 

- ^2 l°° dss ln fyj lPv(s) - pa(s)] , (19) 

where p = 2//e 1//2_c ' = 1.85//. The use of // here is just for convenience; AD is 
actually independent of it. 

3 Phenomenology at Zero Temperature 

Phenomenological analysis of the Weinberg sum rules was originally made in a 
very simple approximation using only the contributions of p,ai,ir mesons. In other 
words, Weinberg assumed that contributions from all excited states other than the 
lowest resonances mentioned cancelled out. Together with the KSFR relation it 
leads to the famous prediction m ai = \f2m p which looked excellent from the point 
of view of data available at the time. However now we know that this prediction, 
as well as predictions for coupling constants, agrees with experiment only up to the 
level of 10-20%. 
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The sum rules are exact in the chiral limit, so one should be willing to verify 
them as accurately as possible. If the complete spectral densities were measured 
one could simply evaluate the integrals and check whether the sum rules are indeed 
satisfied, up to the accuracy of the data. Unfortunately the situation is not that 
straightforward because there are no meaningful measurements of the non-resonance 
contribution in the axial channel. Therefore we first have to close this hole using 
the sum rules themselves. 

Let us first discuss how well the spectral densities are determined experimen- 
tally. In the pole plus continuum approximation one would write 



t s m t r/ 9n S l\ a s( S ) 

py(s) = ^ (5 ( s _ m 2 ) + _ 1 + ^ + 



9(s - E 2 V ) , (20) 



and 

PA(s) = ^5(s - ml) + * [l + ^ + • • •] 9(s - E\) . (21) 

9 a L TT J 

The continuum is the same in the vector and axial-vector channels, according to 
perturbative QCD, but the phenomenological threshhold is in general different. 
Note that the individual integrals over s for the vector and axial-vector channels 
which enter the sum rules are actually divergent because of the continuum, but 
the difference is finite. The coupling constants are the same ones used in a vector 
dominance approximation to the currents as expressed in the current-field identities 
of Sakurai, 

2 

vn 

K = -tpI > ( 22 ) 

9p 

2 

77? 

^ = — < + Pion- (23) 

9a 

We don't use these approximations for the spectral densities because the three sum 
rules involve integrations of the spectral densities with different powers of s and so 
it is likely important to incorporate the finite widths of the resonances. 

The vector spectral density is very well measured in e + e~ — > p — > 7r + 7r~. An 
s-wave relativistic Breit-Wigner is not a good representation because the p meson is 



a p-wave resonance. A much better representation is given by the Gounaris-Sakurai 



formula [17, 18]. It turns out that this complicated formula can be approximated by 
a relativistic Breit-Wigner with an effective width T' p = 118 MeV and an effective 
pole mass m' p = 761 MeV. 



Pv{s) 



4 1 

nip l 



m'p) 2 + nipT'p 2 



+ 



8vr 2 l + exp[(£V- ^fs)/5 v \ 
We take g p from the KSFR relation 



l + 



0.22 



ln(l + yi/0.2GeV) 



4 



m. 



(24) 



(25) 



With nip = 768 MeV and F n = 94.5 MeV one gets from this g^/Atr = 2.63. The 
second term in eq. (24), corresponding to the continuum from 2n-pion states (n = 



2, 3, ...), has E v = 1.3 GeV and by = 0.2 GeV 0. 

The coupling of the a\ to the current can be determined from the measured 
branching ratios of r — > v T + hadrons. According to the 1992 Particle Data Table, 
the branching into the two 3-pion channels dominated by the a\ is ll. 2 ± 1.4% 
while the p dominated 2-pion channel is 24.0 ± 0.6%. The first number gives rise 
to the main uncertainty in our numerical analysis below. Using these numbers and 
the theoretical expression for the branching ratios (which follows from the narrow 
width approximation) 



B(t -^v t + oi) _ m 2 ai g 2 (1 - mjjm 2 ) 2 {l + 2m 2 a Jm 2 ) 



(26) 



B(t — > v T + p) m 2 g 2 ± (l — m 2 /m 2 ) 2 (l + 2m 2 /m 2 ) 
one can get the coupling g ai = 10.5 ± 0.7. 

For the axial-vector spectral density we use an expression analogous to the 
vector one but with the following differences. First, we use a constant width of 400 
MeV and a constant mass of 1260 MeV for the a\ contribution (for more details 
about this see reference fil| ). We have, however, cut off this resonance below the 
threshold m p + m w . 
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The large width of the a\ and its proximity to the r lepton causes a significant 
correction to eq. (26). Numerically integrating the differential decay rate g with 
the realistic shape of the resonance we get finally a value g ai = 9.1 ± 0.7. 

The available data for the nonresonant axial states are very poor so that the 
continuum threshhold Ea and the width 5a are unknown. The reason is partly 
statistical. More importantly, since the data about the axial spectral density come 
from the r lepton decay, there are fundamental limitations due to the r mass which 
is not big enough to provide sufficient phase space for 3- and 5-pion final states with 
the needed invariant mass. Therefore, some authors (for example have analyzed 
the Weinberg sum rules without the axial continuum. 

In Fig. 1 we show our spectral density with axial continuum using the same 
width as the vector continuum [|l]] and with a threshhold value to be determined 
below. In this figure the dashed curves correspond to the experimental uncertainty 
in the branching ratio into a\. One can see that this is a rather non-trivial, sign- 
changing function, which should obey the sum rules under consideration. Naturally, 
the new sum rule we consider is more sensitive to the large s behavior of the differ- 
ence of the spectral densities. Thus we may at least ask whether all sum rules are 
consistent with one common value of the parameter Ea- 

In Fig. 2 (a-c) we have plotted sum rules I to III as functions of Ea- The 
horizontal dashed line shows in all cases the right side of the sum rule which depends 
on the vacuum quark condensate or F v as appropriate. The intersection of the lines 
should occur at the same value of Ea- As explained in the previous section, we do 
not know exactly the VEV of the relevent operators, therefore we use the vacuum 
dominance estimate, with 

| < 0\uu\0 > | 1/3 = 240 MeV (fi = 1 GeV) . (27) 

Fortunately, there seems to be very little sensitivity to the value of the quark con- 
densates. One can clearly see that sum rules II and III are quite consistent with 
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the common value of Ea = 1-45 GeV. This observation is nontrivial. 

Now we can come back to the first sum rule, use this value of Ea as input, and 
compare the numerical value of the integral to the right hand side. This procedure 
predicts F n about 5% higher than the experimental value. 

Finally, let us comment on a closely related integral of the spectral densities 
under consideration. It was shown in |?0| that the electromagnetic mass difference 
of pions can be expressed as 



where A is some cutoff parameter used to regulate the divergent integral over virtual 
momentum in the loop. The result obviously depends on it; only in one particular 
limit, namely, for A 3> m p , m ai and for the original Weinberg values of the p and a± 
parameters without continuum one can get rid of it and recover the original result 



data as explained above, it is no longer true. The integral does depend on the cutoff 
A. 

In Fig. 2 (d) we show this sum rule with A = 2 GeV as a function of Ea- 
Note that the value of the pion mass splitting is very sensitive to Ea, and can even 
change sign if it is only 40 MeV above the suggested value. However, at Ea = 1-45 
GeV it agrees with the experimental value (horizontal line) reasonably well. Fine 
tuning could be accomplished by adjusting the cutoff A, but we shall not do this. 

4 Finite Temperature Sum Rules 

In this section we first generalize Weinberg's two sum rules to finite temperature 
using essentially the same methods as he used without any specific reference to 
QCD. Then we verify the generalizations by using the OPE, which also allows us to 
obtain the finite temperature extension of sum rule III. Finally, we investigate the 




(28) 



vr ^+ — m^lo = (3 In 2a/27r)m 



? p of pC[ |. However, for the parameters extracted from 
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behavior of these sum rules at low temperature. 



4.1 Derivation of Weinberg-type sum rules at fixed momentum 

Consideration of Weinberg-type sum rules at finite temperature (or chemical poten- 
tial) is more involved than at zero temperature. Lorentz invariance is not manifest 
because there is a preferred frame of reference, the frame in which the matter is at 
rest. Thus spectral densities and other functions may depend on energy and mo- 
mentum separately and not just on their invariant s. Also, the number of Lorentz 
tensors is greater because there is a new vector available, namely, the vector = 
(1,0,0,0) which specifies the rest frame of the matter. 

For a given 4-momentum p it is useful to define two projection tensors. The 
first one, P£" , is both 3- and 4-dimensionally transverse, 

pn = 5 ij _ t?L ; (29) 

P 

with all other components zero. The second one, P^", is only 4-dimensionally 
transverse, 

P£ v = - (V - + Pf) ■ ( 3 °) 

The notation is L for longitudinal and T for transverse with respect to p. There 
are no other symmetric second rank tensors which are 4-dimensionally transverse. 

We now define the longitudinal and transverse spectral densities for the vector 
current as 

< V£(z)V?{0) >=J^J d 4 P 9(p°) \&P£ + p^Pf] , (31) 
and for the axial vector current as 

< i4£(*K(0) >=J^sj d 4 P e(p°) e** [p L A Pr + p T A Pf] . (32) 

In these expressions the angular brackets refer to the thermal average. In general 
the spectral densities depend on p° and p separately as well as on the temper- 
ature (and chemical potential). These definitions are standard and insure that 
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both spectral densities are non-negative. In the vacuum we can always go to the 
rest frame of a massive particle, and in that frame there can be no difference 
between longitudinal and transverse polarizations, so that Pl = PT = P- Since 
P^ + = -(g» v - p^p u /p 2 ) these equations collapse to eqs. (8-9). The pion, 
being a massless Goldstone boson, is special. It contributes to the longitudinal axial 
spectral density and not to the transverse one. In fact, we could write 



(33) 



This shouldn't be done at finite temperature because the contribution of the pion 
to the longitudinal spectral density cannot be assumed to be a delta function in 
p 2 . In general the pion's dispersion relation will be more complicated and will 
develop a width at nonzero momentum. Therefore, we do not try to separate out 
the pionic contribution but subsume it in the spectral density p^, without any loss 
of generality. 

Following Weinberg, we define a three-point function by 



ie abc M^\q,p) = J d 4 xd 4 ye 



-i(q-x+p-y) 



A»{x)At{y)V c \Q) 



> 



(34) 



We multiply both sides with q^. On the right side we can use 



-i(q-x+p-y) 



, d 



i{q-x+p-y) 



(35) 



Both the vector and axial-vector currents are conserved. We assume that we can 
integrate by parts and that the surface term is zero. The nonzero contribution 
comes from 

^{T [A^x)Al{y)V c \Q)]} 
= 5(x -y°){9(x ) [Al{x),Al{y)] V c x (0) + e(-x°)V c x (0) [^(x), j4*(y)] } 



+ 



5(x°) {e(y°)AUy) [A° a (x),V c \0)] +9(-y°) [A° a (x), V c x (0)] A» b (y)} . (36) 
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From this expression we see the need for knowledge of the equal time commutators. 
Consistent with the normalization of eqs. (5-7) we have 



5(z°)\A° a (x),AUy) 



d 

ie abd V£(x)8(z) + Sy J ab (x)-^-j5(z) 



5(z°)\A a (x),V b »(y)] = ie abd A^x)5(z) + SZ b (x)g-5(z) 



(37) 



Here z = x — y, and the <S"s denote the Schwinger terms. 

Consider now the contribution of the Schwinger terms to the thermal average. 
Generically they will be of the form 



< SJ>= Z- l Y,e 

m,n 



-K n /T 



< n\S\m >< m\J\n > , 



(38) 



where K = H—fiN is the Hamiltonian minus the chemical potential times conserved 
particle number, the states are chosen to be eigenstates of H, N, and isospin, and 
J is either the vector or the axial-vector current. J has isospin one, so we get zero 
if either (i) S is a c-number, or (ii) S is an operator with no isospin one component. 
We shall assume that one of these holds. Then 

d 



dx^ 



> 



ie abd S(x -y)<T [V7(x)V c A (0)j > + ie acd 6(x) < T [A^y)A^0) 
It is now a simple matter to show that 

l -q»M^\q,p) = Dp(q+p)-Df(p), 

where the D's are the propagators for the currents, such as 



> 



(39) 



(40) 



\p) = J d 4 ye' ip y <T 
Similarly, one can show that 



<(y)A*(0) 



> 



(41) 



l(q+p)xM^ X (q,p) = D%{q)-D%{p). 



(42) 
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These Ward identities have exactly the same form as at zero temperature [ffl. 
With a similar consideration of the three-point function 

- ie abc N^ x (Q,P) = f d 4 xd 4 ye-^ x+ ^ < T \V£(x)Vf{y)V c x (oj\ > , (43) 



(44) 



one can prove two more Ward identities, 



1 



-q,N^ X (q,p) = Dp(q+p)-Dp(p), 



and 



1 



(q+p) X N^ x (q,p) = D^(q)-D^(p). 



(45) 



Multiply eq. (42) by (q + p)\ and eq. (44) by q^. Do the same for the other 
two Ward identities. One obtains the constraints 



{q+p)xDp{q+p) = qxD^(q) + p x D^(p) = q X D A \q) + p X D v A \p) . (46) 



This implies linearity in the momentum, 

k\Dy X (k) = k x Df(k) = C vX k x , 



(47) 



where C uX is momentum-independent (but can depend on temperature) and is the 
same for the vector and axial-vector channels. By taking the Fourier transform of 
these relations we can find the thermal average of the equal time commutators, 

d 



5(x°)< V:(x),V b °(0) >=5(x°)< A» a (x),A° b (0) 



>= SabC 



i/A 



dx x 



6{x) . (48) 



The commutators above can be expressed in terms of the spectral densities from 
eqs. (31-32). Taking their difference one obtains the finite temperature generaliza- 
tion of the first Weinberg sum rule, 



I 



o to 2 - p 2 



p v (u),p) -p A (u,p) 



0. 



(49) 



Notice that this sum rule involves only the longitudinal spectral densities and not 
the transverse ones. At zero temperature the spectral densities depend only on 
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p 2 = s = uj 2 — p 2 . Then this equation reduces to eq. (1) once we remember to 
separate out the pion piece of p\, namely, sF 2 5(s). At finite temperature, the 
spectral densities in general will depend on uj and p separately and not just on the 
combination s. Then this sum rule must be satisfied at each value of the momentum. 

At this point, Weinberg made an additional assumption in order to obtain the 
second sum rule (eq. (2)): the currents behave like free fields as p 2 — ► oo. He also 
related the difference between the vector and axial-vector propagators to the matrix 
element of a particular operator between the vacuum and a one pion state. This 
is difficult to generalize to an ensemble average. To obtain the finite temperature 
generalization of the second sum rule we follow the arguments of Das, Mathur and 
Okubo pi]] instead. 

Deleting the index V or A the explicit expressions for the propagator and the 
Schwinger term are 




(50) 




(51) 



(52) 



where 




(53) 



(54) 



(55) 



and 



' 00 = C i 



C j0 = 



(56) 



where 




(57) 
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The first observation we can make concerns the thermally averaged Schwinger 
term C. Since it is the same for the vector and the axial-vector correlators, by eq. 
(47), the -Ds(p) must be the same as well. Equating them reproduces the first finite 
temperature sum rule (eq. (49)). 

The essence of the argument of Das, Mathur and Okubo is that spontaneous 
chiral symmetry breaking is a low energy phenomenon. At very high energy it must 
disappear, at least in the limit that quark masses are zero and chiral symmetry is 
exact. Thus the difference between the vector and axial- vector propagators should 
go to zero at very high energy, 



lim 

p°— »oo, p fixed 



D^(p°,p) - D%(p°,p) = 0. (58) 



If we do this for the time-time or time-space components of the propagators, that is, 
for the Dl, we again reproduce the first finite temperature sum rule. Expanding to 
the next order in \jp\ we obtain a finite temperature generalization of the second 
zero temperature sum rule, which is, 

II-L / duu Pv(lo,p) - Pa(^,p) =0. (59) 
Jo 1 J 

Like the first, this sum rule involves only the longitudinal spectral densities, and we 
call it II-L. Also like the first, it reduces to the original Weinberg sum rule as the 
temperature and/or chemical potential go to zero. 

Next we consider the space-space components of the propagators. Examination 
of the D' L in the infinite energy limit gives us the sum rule II-L and nothing new. 
Examination of the Dt in the infinite energy limit gives us another sum rule which 
we call II-T because it involves the transverse spectral densities, 

II-T / duu Pv(lo,p) - Pa{uj,p) =0. (60) 
Jo 1 J 

The finite temperature sum rules II-L and II-T should become degenerate at p = 
because there ought not to be any difference between longitudinal and transverse 
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excitations at rest. The sum rule II-T also then reduces to the original second sum 
rule in the vacuum. 

We want to emphasize that the sum rules derived in this section, I, II-L and II- 
T, must be satisfied for every value of the momentum. Furthermore, our derivation 
is more general than QCD; any theory which satisfies the assumptions we made 
must obey these sum rules. Perhaps they would be useful in the context of models 
of the electroweak interactions where the Higgs particle is a composite of other fields 
or for technicolor theories. 

4.2 Sum rules and the operator product expansion 

Application of the OPE to finite temperature has a peculiar history. In the first 
papers (|| and several later ones) the authors considered only the T-dependence of 
average values of the same operators as at T = 0, the Lorentz scalars. However, the 
rest frame of the heat bath selects a 4-vector, thus symmetric tensors should also 
be included. In fact, the situation is completely analogous to that in deep-inelastic 
scattering, for which one also has a preferred frame, that of the target. Thus, one 
can simply use formulae derived in that context (see discussion in The finite 

temperature sum rules were recently re-examined along these lines in 

The fact that we are not going to discuss vector and axial channels as such, 
but only concentrate on their difference, brings in significant simplifications. Most 
operators describing the interaction of a quark with the gluonic field are chirality 
blind and therefore cancel. In the chiral limit, the difference appears only starting 
with the four-quark operators. 

To leading order in the momentum the difference between the vector and axial- 
vector correlators is given by the OPE to be 




(61) 
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where the operator O was defined in section 2. This structure first appeared in 



the OPE analysis of the next-twist correction to deep inelastic scattering in [22]. 
Observe that this quantity is transverse: p^AD^ v = p u AD^ u = 0. This is con- 
sistent with eq. (47), the equality of the Schwinger terms, and therefore with the 
assumptions made to derive it. 

First, consider AD 00 . In terms of the spectral densities it is given by eqs. (50) 
and (53). Expand it in inverse powers of p\ in the limit that \p§\ — * oo. Since the 
coefficients of 1/pq and 1/pf, in eq. (61) are zero it must be that the corresponding 
coefficients in eq. (50) are also zero. This gives us the finite temperature sum rules I 
and II-L immediately. We can say nothing about the next term without knowledge 
of higher dimension operators in the OPE, which would contribute to order 1/pq. 

Next, consider AD|J. From eqs. (50-55) it is 

duiu 
/o P5 - uj 2 



AD£ = 2i r 2 d0J z . [2Ap T ( W , p) + Ap L (u;, p)l . (62) 



Again, expand in inverse powers of Pq. The term of order 1/pq, when combined 

with the just derived sum rule II-L, gives us the sum rule II-T. The term of order 

1/Pq gives us the finite temperature version of sum rule III. 

f°° r i r i 

III / diou 3 2Ap T (u;,p) + Ap L (uj,p) = -27TQ, < > +2 < O 00 > . 

(63) 

We can make two observations about this sum rule. In the limit of vanishing 
temperature, Lorentz covariance says that 

<0 00 >t=o= \<0»> T=Q . (64) 

This reduces eq.(63) to the previously derived zero temperature sum rule eq. (16). 
At finite temperature, the right side of eq. (63) depends on T but not on p. 
Therefore, the integral on the left side must be momentum- independent. If the 
integral is known at zero momentum, for example, then it must have the same 
value for any momentum. 
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4.3 The low temperature limit 

As we are taking the zero quark mass limit in this work, the pion is massless below 
any critical temperature for chiral symmetry restoration and/or deconfinement, and 
thus at parametrically low temperature the heat bath is dominated by pions. In 
[p3| the so-called Dey-Eletsky-Ioffe mixing theorem was proven, which says that, to 
order T 2 , there is no change in the masses of vector and axial- vector mesons. What 
changes are the couplings to the currents. The finite temperature correlators can be 
described by a mixing between the vector and axial-vector T = correlators with 
a temperature dependent coefficient, 

D^(p,T) = (l-e)D^(p,0) + eZ>7(p,0), (65) 
D%(p,T) = (l-e)D^(p,0) +eD(T(p,0). (66) 

These are valid to first order in e = T 2 /6F%. This implies the same mixing of the 
spectral densities, namely, 

PV (P°,P,T) = (l-e) Pv (s,0) + ep A (s,0), (67) 
p A (p°,P,T) = (l-e) PA (s,0) + ep v (s,0), (68) 

with the appropriate longitudinal and transverse subscripts. The temperature de- 
pendence of the pion decay coupling was thus proven to be F%(T) = (1 — e)F% for 
small T consistent with the prediction of chiral perturbation theory j24|. Therefore, 
the finite temperature sum rules I (eq. (49)), II-L (eq. (59)) and II-T (eq. (60)) 
reduce to the original, zero temperature sum rules but with both sides of the eqs. 
(1) and (2) multiplied by the factor 1 — 2e. 

One may ask whether the third sum rule also obeys the Dey-Eletsky-Ioffe mixing 
theorem. A general formula describing the thermal average of any four-quark oper- 



ator using soft pion methods was derived in [25]. For an operator Oab = qAqqBq 
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the expression is 



<Oab > 



144 ^ 

< uu > 2 T 2 
144 12F2 



< uu > 2 




4F 2 



T 2 



) 



[Tt(A)Tt(B) - Tr(AB)} 



[Tr( l5 T a A)Tv( 75 T a B) - Tx( l5 r a A l5 T a B)} , (69) 



where it is assumed that at T = one can use the vacuum dominance approximation. 

The average of the four-quark operator appearing in sum rule III gets multiplied 
by the correct factor 1 — 2e, as shown by Eletsky [p5f| . This is not a trivial result: 
the average value of an arbitrarily chosen four-quark operator will not have the 
same temperature dependence. As already emphasized by Eletsky, a simplistic 
application of factorization at nonzero temperature, which would suggest the same 
behavior as for the quark-condensate squared, 



would be wrong, and in fact violates the sum rule. 

In summary, at low temperature the sum rules under discussion satisfy the 
Dey-Eletsky-Ioffe mixing theorem exactly. 

5 Scenarios for Chiral Symmetry Restoration 

Chiral transformations are rotations of the quark field with 75 , and they may or may 
not have the SU(Nf) (isospin) generators. The corresponding U(1)a and SU(Nj)a 
have different fates in QCD; the former is explicitly violated by the anomaly, the 
latter is broken spontaneously at low temperature and is restored at some critical 
temperature T c , provided the quark mass is stricly zero as it is assumed in this 
paper. However, as the p and a\ channels are both isospin-1, the symmetry which 
can mix them is U(1)a- 

If both chiral symmetries are restored, one may conclude that left and right- 
polarized quarks form completely independent subsystems. If so, quarks never 




(70) 
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change chirality, and there should be no difference between vector and axial-vector 
correlators, and the object of our consideration is simply zero. 

This is indeed expected to happen at very high temperatures, but in the critical 
region T ~ T c all depends on the mechanism of the U(1)a symmetry breaking. One 
of us recently wrote a mini-review on this subject p2| , and it is probably enough to 
mention here that it is most likely that the U{\)a symmetry is practically restored 
at T ~ T c . More specifically, there is direct evidence from lattice simulations that 
the difference between vector and axial-vector correlators do indeed vanish around 
this point, to within the accuracy of the calculation. 

In this section we speculute on exactly how this difference goes to zero with 
increasing temperature. Generally, one may suggest many different scenarios. Let 
us discuss the following three. 

5.1 Mixing of vector and axial- vector spectral densities 

The simplest scenario is that the T-dependence factorizes. It means that the vector 
and axial-vector spectral densities mix, without changing their shape, as in the 
low temperature limit considered in the previous section, only with a more general 
function e(T). When the mixing becomes maximal, e = 1/2, chiral symmetry is 
restored. 

It is amusing to see at what temperature this occurs using the lowest order 
formula, e = T 2 /6-F 2 . This estimate gives T com p\ e te mixing = V3F n « 164 MeV, 
which is indeed roughly equal to the expected critical temperature T c . 

5.2 Shift in meson pole masses and residues 

In this scenario we assume that the p and a\ mesons retain their identities and dom- 
inate the correlation function. However, their parameters change with temperature. 
In particular, the masses may move towards each other ||, or go to zero ||. At T c 
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they become degenerate, and chiral symmetry is restored. 

It is instructive then to look at the sum rules. Let us assume that vector meson 
dominance is a good approximation for the spectral densities and not worry about 
the continuum contribution for the time being. Let us focus on zero momentum 
for the sake of simplicity. When a pole mass is defined at finite temperature, it is 
usually defined as the energy of the excitation at zero momentum. 

The vector spectral density is (there is no difference between longitudinal and 
transverse at zero momentum) 

1 TTl^ 1 

p v (uj) = ^hm — p p , (71) 

vr g 2 u 2 -m 2 - U r (uj) - iTVj{uj) 

where 11^ and Ii p j are the real and imaginary parts of the p self-energy at tempera- 
ture T. In the narrow width approximation this becomes 

pv{u) = His^-ml-U^u)) . (72) 

The pole mass is determined self-consistently from m 2 (T) = m 2 p +Jl p R {rn p {T)). Then 
the spectral density can be rewritten as 

p v {u) = Zp(Ty^-5 (u? - m 2 p(T)\ , (73) 



where the temperature-dependent residue is 

d 



Zp-\T) = 

The normalization is Z p (0) = 1. Similarly 



1 " ^ n «<"" 



(74) 



4 

p A (v) = Z a {T)^5 (u, 2 - m 2 ai (T)) + Z w (T)F 2 uj 2 5 (u 2 ) . (75) 

Substituting these spectral densities into the finite temperature sum rules I and 
II-L/II-T tells us that the p and a\ residues are equal 

Z p {T) = Z a {T) , (76) 
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and that the pion residue is 



Z n (T) = 2Z p {T) 



2 1 



(77) 



mj(T) ml(T)_ 

We expect that m^(T) — m^(T) — > as the temperature increases. Three types of 
behavior can be distinquished: both the p and the a\ masses decrease with T, both 
masses increase with T, or the p mass increases while the a\ mass decreases with 
T. The sum rules do not appear to rule out any of these possibilities. In any case, 
the result is that Z n (T) — > unless Z p (T) — > oo, which seems rather unphysical. 

5.3 Resonance broadening and downward shift of the continuum 

As distinct from the previous scenarios, it may be that particles are not well-defined 
as we approach a chiral symmetry restoring phase transition. That is, the imagi- 
nary part of the self-energy may become larger with increasing temperature. This 
broadening would also decrease the maximum peak value of the spectral density. 
Euphemistically, the vector and axial-vector mesons melt away. There may be also 
a decrease in the threshholds Ey{T) and Ea{T) of the continuum. The continuum 
would merge with the broadened particle poles to give a very broad distribution 
of strength in the spectral densities. The difference of spectral densities shown 
for T = in Fig. 1 would become flatter and decrease everywhere towards zero, 
effectively restoring chiral symmetry. 

Concluding this section, we say once more that the sum rules by themselves 
cannot of course tell which scenario is preferable. However, the sum rules can be 
used to significantly restrict the parametrization of the spectral densities at nonzero 
temperature. 
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6 Conclusion 



In this paper we studied Weinberg-type sum rules at zero and at nonzero tempera- 
ture. All considerations were made in the exact chiral limit of QCD, m q — > 0. In the 
former case, we derived a new sum rule of the Weinberg-type. Although it belongs 
to an infinite series of sum rules, one for each type of OPE term at small distances, 
we think it is special in several respects. First, it is relatively simple theoretically 
because it is related to the VEV of a specific four-quark operator. Sum rules of 
higher order than the third are much more complicated. Second, it is related to 
the leading nonzero ln(r) term of the correlators, while others can be related to 
sub-leading terms which are much more difficult to single out, especially in lattice 
simulations. 

Continuing the zero temperature analysis, we re-examined the experimental 
data together with all relevant sum rules. We found that, although we do not have 
sufficient information on the VEV of the operator for sum rule III, we still can use it, 
together with sum rule II, to fix the numerical value of Ea, the continuum threshold 
in the axial channel. This essentially closes the gap in the experimental data, and 
allows one to test the original Weinberg sum rules without any ad hoc assumptions. 
Good agreement with the experimental values of F w and the electromagnetic mass 
difference of pions provides a non-trivial consistency check of the data used. 

Our finite temperature analysis consists of several different parts. First, follow- 
ing Weinberg's original derivation, one can find generalizations of his sum rules to 
nonzero temperature. Sum rule I involves only the difference of the longitudinal 
spectral densities, while sum rule II bifurcates into two sum rules, one involving the 
longitudinal spectral densities, and the other involving the transverse ones. These 
sum rules must be satisfied at each value of the momentum. These new features 
arise because of the appearance of a preferred reference frame at nonzero tempera- 
ture. These sum rules were derived without specific reference to QCD so they are 
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applicable to other theories satisfying the assumptions made. We were also able 
to derive them from the OPE. Furthermore, we used the OPE to obtain the finite 
temperature generalization of the new sum rule III, which makes specific reference 
to the dynamics of QCD. 

We also considered very low temperatures at which chiral perturbation methods 
predict the general behavior of the correlators. We showed that these results are in 
exact agreement with all sum rules under consideration. 

We would like to emphasize that the average value of the four-quark operator 
which appears in sum rule III is of great theoretical interest. It shows correlation 
between densities and currents made of left and right-handed quarks and is, in a 
sense, an order parameter for restoration of U (I) a chiral symmetry. The average 
value in the QCD vacuum and and at finite temperature can and should be studied 
in lattice numerical simulations. This task is facilitated by the fact that it does not 
have any perturbative contributions. 

Finally, we speculated on possible scenarios of chiral symmetry restoration. We 
have no preferences among them, and only future work, including especially lattice 
numerical simulations, can clarify which of them (if any) is realized in QCD. How- 
ever, the derived sum rules should hold in any case, thus providing some relations 
among parameters of the vector and axial-vector spectral densities. 
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Figure Captions 

Fig. 1. The difference between the spectral densities of the vector and axial-vector 
currents versus s at zero temperature. The two dashed curves show the uncertainty 
due to the experimental determination of the a\ coupling constant, described in the 
text. The abrupt change at s = 0.8 GeV 2 corresponds to the sharp onset of the a± 
contribution at [m p + m n ) 2 . 

Fig. 2. Dependence of the zero temperature sum rules I-III on the effective pertur- 
bative threshhold Ea in the axial channel. The last panel shows the ir + — ir° mass 
difference. As in Fig. 1, the solid and the two dashed curves correspond to the 
central value and uncertainty in the a\ coupling constant. In all cases, the expected 
magnitude of the corresponding sum is shown by the horizontal long-dashed line. 
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